30th International Cosmic Ray Conference 




ICRC'^/ 



Merida, Mexico 



1^ 

o 
o 

(N 



Photon spectra from final stages of a primordial black hole evaporation in different 
theoretical models 

Edgar Bugaev, Peter Klimai, Valery Petkov 
Institute for Nuclear Research, Moscow, Russia 
bugaev @ pcbailO. inr. ruhep. ru 

Abstract: Possibilities of an experimental search for gamma-ray bursts from primordial black hole (PBH) 
evaporations in space are reconsidered. It is argued that the corresponding constraints which can be 
obtained in experiments with cosmic ray detectors strongly depend on theoretical approach used for a 
description of the PBH evaporation process. Predictions of several theoretical models for gamma-ray 
spectra from final stages of PBH life (integrated over time) are given. 
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Introduction 

Direct searches for the bursts of gamma rays from 
the evaporations of PBHs have been carried out in 
several works during last decade [2, 3, 6, 8, 10]. 
PBHs are a unique probe of general relativity, cos- 
mology and quantum gravity [4]. They could have 
been formed in the early universe, e.g., from di- 
rect gravitational collapse of primordial density 
perturbations. It was shown in numerous works 
that PBHs can be produced due to specific features 
of the inflationary potential or due to perturbation 
amplification during preheating phase. For direct 
searches of PBHs their spatial clustering proper- 
ties are very essential. It was shown recently [5] 
that the local clustering enhancement of PBH num- 
ber density in our galaxy can be very large, up 
to ~ 10^^, i.e., many orders of magnitude larger 
than the factors ^ 10^ predicted earlier in some 
works. Due to this, the limits from direct searches 
of PBHs can be much stronger than indirect lim- 
its from measurements of diffuse gamma ray ex- 
tragalactic background. 

PBH evaporation models 

Model of MacGibbon and Webber 

In first calculations of photon spectra from BH 
evaporations it had been shown that this spectrum 



is far from being thermal simply because emitted 
elementary particles such as quarks fragment into 
hadrons, photons, neutrinos, etc, and just this frag- 
mentation and decay of unstable hadrons produces 
the final photon spectrum [11]. It was assumed that 
the radiation evaporated from the BH interacts too 
weakly and all emitted quarks propagate freely and 
fragment independently of each other 

The resulting time-integrated spectrum (for the 
case Th ^ m-^) can be parametrized in the simple 
form 
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where Eq — 10^, and all energies are measured 
in GeVs. The parametrization is valid for E^ > 
m7r/2. The average energy of this distribution is 
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Ej « 2Th- 
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« 1300At~i/^ , (1) 

where we used, just for example, Emin — 
lOOGeV, and At is the remaining PBH lifetime (in 
seconds), which is connected with initial Hawking 
temperature by the simple relation: 



At ^ 4.7 X 10' 
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Figure 1: Instantaneous photon spectra from a 
black hole with Hawking temperature Th = 10 
TeV. 
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Figure 2: Time-integrated photon spectra from a 
black hole with initial Hawking temperature Th = 
10 TeV. 



Heckler model 

According to Heckler's idea [9], once a black hole, 
in course of an evaporation, reaches some crit- 
ical temperature Tcru, the emitted Hawking ra- 
diation begins to interact with itself and forms a 



nearly thermal chromosphere. This assumption is 
based on a simple argument: at relativistic ener- 
gies, the cross section for gluon bremsstrahlung 
in a qq or qq colUsion is approximately constant 
(up to logarithms of energy) while the density of 
emitted particles around the BH is proportional 
to BH temperature. The evaporation of BH cre- 
ates a radiation shell which propagates outward 
at the speed of light. The quark and gluon spec- 
trum in the observer rest frame is obtained by 
boosting a thermal spectrum at Tcru = ^qcd 
(and Aqc d = 200McV) with the Lorentz factor 
7cft » 0.22{Th/Aqcd)^^'^ of the outer surface of 
the thermal chromosphere. 
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In this formula the integration over the surface of 
the chromospere of radius rch = ^chj Aqcd has 
been carried out. The observed photon spectrum 
of the chromosphere is a convolution of the quark- 
gluon spectrum given by (3) with the quark-pion 
fragmentation function and the Lorentz-boosted 
spectrum from tt" decay. 

The high energy behavior of time-integrated spec- 
trum is {E^ is in GeVs) 



(4) 



The average energy of time-integrated spectrum is 
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Kapusta & Daghigh model 



This approach [7] differs from those of Heckler 
mainly in two respects. It is assumed that the 
hadronization of quarks occurs before the thermal 
freeze-out and the onset of free- streaming. It is as- 
sumed, further, that it happens suddenly at a tem- 
peratute T/ in the range 100 - 140 MeV (it is the 
parameter of the model), and to this moment all 
particles whose mass is greater than T/ have an- 
nihilated leaving only photons, electrons, muons 
and pions. Photons emitted by BH come from two 



30th International Cosmic Ray Conference 



sources. Either they are emitted directly from a 
(boosted) black-body spectrum or they arise from 
the 7r°-decay. 

The formula for particle spectrum from the outer 
edge of the chromosphere is analogous to (3). For 
the concrete calculation, we used Tf = 0.12 GeV. 
The values of jch and rch obtained in [7] are 
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The gamma ray flux from the 7r° - decay in the limit 
^ rTijr is given by the simple formula [7] 



(7) 



and we used it to estimate photon flux from this 
channel even at ^ m.^. 

In this case, the high-energy behavior of time- 
integrated photon flux is similar to (4), but the flux 
is approximately 1 .5 orders of magnitude larger: 
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The average energy is given by 

0.36VTh" « 32At"i/^ 
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Experimental search of PBH bursts 
with ground based detectors 
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Figure 3: Average energy of gamma rays emitted 
from a PBH, as a function of left PBH lifetime. 
Upper curve corresponds to Emin = lOOGeV. 
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Figure 4: Average photon energy vs. zenith angle 
for model of [7] . 



One can see from Fig. 3 that the average pho- 
ton energy as a function of PBH lifetime is highly 
model dependent. Therefore, various methods 
should be used for the search. Up to now the search 
of high energy gamma radiation from PBH with 
ground-based arrays had been carried out using the 
predictions of the model of [11], and, specifically, 
data on distribution of events (extensive air show- 
ers (EAS)) in time and space were used. The re- 
sults were reported in [2, 3, 6, 8, 10], time intervals 
used were from 1 to 10 seconds. 

If the chromosphere models of black hole evapo- 
ration are close to reality the average photon en- 
ergy at last instants of PBH life is rather low. In 
the model of hadronic chromosphere developed in 
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Figure 5: Time duration of gamma burst vs. zenith 
angle for model of [7]. 



Photon spectra from final stages of a primordial black hole evaporation 



[7] this energy is about 30 GeV for At ~ Is. A 
search of PBH bursts in this case is possible using 
the method suggested in works of EAS-TOP col- 
laboration [1] and Baksan group [12, 13]. Shower 
particles generated in atmosphere by photons with 
energy ~ 30 100 GeV are strongly absorbed 
before reaching the detector level, so, the average 
number of signals in the detector module is smaller 
than 1 . Therefore, in this energy range PBH bursts 
can be sought for by operating with the modules 
in single particle mode, that is, by measuring the 
single particle counting rate of the individual mod- 
ules. The primary arrival directions of photons are 
not measured, and bursts can be detected only as 
transients (short-time increases) of the cosmic ray 
counting rate. 

Using the photon spectra from PBH bursts calcu- 
lated above and registration efficiencies of photons 
(as functions of zenith angle, the array altitude and 
features of concrete modules) one can calculate the 
average energies of PBH bursts for a given array 
and, approximately, corresponding time durations 
of these bursts. 

On Fig. 4 the dependencies of average photon en- 
ergy on zenith angle in model of [7] for two arrays 
of Baksan observatory are shown. Andyrchy and 
Carpet arrays are situated practically on the same 
place (horizontal distance between them is about 1 
km), but their altitudes above sea level are differ- 
ent: Andyrchy - 2060 m. Carpet - 1700 m. On Fig. 
5, time duration of bursts as a function of zenith an- 
gle is given for the same arrays. One can see that 
duration of burst strongly depends on the zenith an- 
gle. 

In summary, it is shown that predictions for high 
energy photon spectra from PBH burst depend 
rather strongly on the model used for description 
of final stage of PBH evaporation. In particular, 
chromosphere models predict too steep spectra in 
TeV region. It is argued that, if such models are 
correct, the practical search of PBH bursts can be 
carried out using the technique developed by EAS- 
TOP and Baksan groups for a search for 7-bursts 
at energies -B > 10 GeV. It is shown that in such 
search the average energy and time duration of the 
signal depend on the zenith angle of its arrival and 
on the altitude of the array above sea level. 

The work was supported by Russian Foundation 
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